Experimental measurement of acceleration correlations and pressure structure 
functions in high Reynolds number turbulence 
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We present measurements of fluid particle accelerations in turbulent water flows between counter- 
rotating disks using three-dimensional Lagrangian particle tracking. By simultaneously following 
multiple particles with sub-Kolmogorov-time-scale temporal resolution, we measured the spatial 
correlation of fluid particle acceleration at Taylor microscale Reynolds numbers between 200 and 
690. We also obtained indirect, non-intrusive measurements of the Eulerian pressure structure 
functions by integrating the acceleration correlations. Our experimental data provide strong support 
to the theoretical predictions of the acceleration correlations and the pressure structure function in 
isotropic high Reynolds number turbulence by Obukhov and Yaglom in 1951. The measured pressure 
structure functions display K41 scaling in the inertial range. 
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Fluid particle acceleration is an important quantity in 
turbulent flows For example, it plays a significant role 
in the formation of cloud droplets in the atmosphere . 
In recent years, advances in the study of the statistics of 
acceleration have been made through the development of 
Lagrangian experimental techniques [1, 0, Q and the use 
of direct numerical simulations (see, e.g., Refs. [1, 
Using silicon-strip detectors operating at recordingfre- 
quencies as high as 70 kHz, La Porta and co-workers [3, Q 
were able to follow passive tracer particles in a water flow 
at Taylor microscale Reynolds numbers up to Rx ^ 10^. 
Fluid particle accelerations were obtained from the tra- 
jectories of the tracer particles. Their results revealed 
the highly intermittent nature of acceleration, and also 
showed the necessity of sub-r^ temporal resolution for ob- 
taining accurate acceleration measurements, where is 
the Kolmogorov time scale, the smallest time scale in tur- 
bulence. In later studies, the same technique was used to 
investigate the Lagrangian properties of acceleration fol- 
lowing a fluid particle [g, [§] . Due to the one-dimensional 
nature of the silicon-strip detectors, however, only one 
fluid particle could be followed at a time. Consequently, 
the spatial properties of acceleration were not explored 
in these previous studies. In other particle tracking ex- 
periments, digital cameras were used to record the mo- 
tion of tracer particles and multi-particle statistics were 



obtained jld.lll|. These experiments, however, were lim- 



ited to flows with small Reynolds numbers because of 
the slower recording frequency of the cameras. Very re- 
cently, advances in CMOS camera technology have pro- 
vided the opportunity of measuring the acceleration of 
multiple tracer particles simultaneously in high Reynolds 
number turbulent flows [l^ . In this Letter, we present 
the flrst direct experimental measurement of the spa- 
tial correlations of acceleration in turbulent flows with 



200 < i?A < 690. 

Another important quantity in high Reynolds number 
turbulence that is not clearly understood is pressure. It 
has been shown that the clustering of inertial particles in 
turbulence is related to the scaling properties of the pres- 
sure field [13] ■ It is, however, extremely difficult to mea- 
sure pressure in turbulent flows non-intrusively. Ould- 
Rouis et al. [11] reported that, in the inertial range, the 
pressure structure functions computed from the fourth 
order longitudinal velocity structure functions scale as 
predicted by Kolmogorov's K41 theory [3, 13, [l3] when 
the Reynolds number is moderately high {R\ > 230). 
However, Hill & Boratav [l^ argued that very large 
Reynolds numbers are needed to observe K41 scaling and 
the assumptions made by Ould-Rouis et al. result in large 
uncertainties in the calculated pressure structure func- 
tion. Pressure spectra obtained from numerical simula- 
tions [13] suggested that the K41 pressure spectrum can 
only be observed at R\ > 600. The spectra obtained 
from direct pressure measurements in turbulent jets by 
Tsuji & Ishihara [20| seem to support this conclusion. 
In this experiment, however, the effect of Taylor's frozen 
flow hypothesis on the pressure spectra has not been fully 
evaluated. 

In high Reynolds number turbulence, the acceleration 
is mostly determined by the pressure gradient, and the 
viscous term may be ignored 2l| . Under this assumption, 
there exist analytical relations between the spatial cor- 
relations of acceleration and the Eulerian pressure struc- 
ture function 3, 17]. By exploiting such relations, we 
obtain an indirect but non-intrusive measurement of the 
pressure structure functions in high Reynolds number 
turbulence. Our experimental results strongly support 
the theoretical predictions of Obukhov & Yaglom based 
on K41. 
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In homogeneous, isotropic turbulence, the Eulerian 
pressure structure function depends only on the sepa- 
ration distance, i.e., n(r) = {[p{x + r) — p{x)]'^) = n(r). 
For high Reynolds number flows, the fluid acceleration is 
dominated by the local pressure gradient. Hence, 



1 / dp 

i?y(r) EE (a,(x)aj(x + r)) = ^(^g^ 



dp 



x+r 



(1) 



Using the relation between the correlation of pressure 
gradients and the pressure structure function in homo- 
geneous turbulence, the acceleration correlation can be 
related to the pressure structure function as [16j 



%(r) 



1 <92n(r) 



2/3^ dridrj 

In homogeneous, isotropic turbulence, this reduces to 
1 d2n(r) 



(2) 



RLL{r) 



2p2 dr2 



RNN{r) 



1 dn(r) 



(3) 



where RLL{r) and RNNir) are the longitudinal and 
transverse acceleration correlations, respectively. There- 
fore, once either n(r) or Rij{v) is determined, the other 
can also be obtained. It should be emphasized that 
Eqs. ([3]) hold in homogeneous, isotropic turbulence at 
high Reynolds numbers. The only simplification invoked 
in deriving these equations is the neglect of a viscous 
contribution to acceleration. 

In their work, Obukhov & Yaglom further assumed, 
as first proposed by Millionshchikov [22| and Heisen- 
berg [2^, that the components of the velocity gradient 
are drawn from a multi-dimensional Gaussian distribu- 
tion. Under this hypothesis, the pressure structure func- 
tion in homogeneous, isotropic turbulence satisfies 



dr4 



(4) 



where $(r) can be written in terms of the derivatives of 
the longitudinal velocity structure function DLL{r): 



6D'r 



LL 



(5) 



An equation for D^i can be obtained from the Karman- 
Howarth equation as 



+ \S\[DLL{r)\ ^ Tsr, 



dr 



(6) 



where S is the structure function skewness [25| . Obukhov 
& Yaglom [3] assumed that S is constant for all sepa- 
rations r, and so it can be related to the Kolmogorov 
constant C2 for the structure function D]^i^{r) as |5| — 



-3/2 

from experiments [26 



(4/5)C2 '''^ and the value of C2 = 2.13 is well known 
Upon solving Eq. ^ numerically 
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1: (Color online) Measured normalized acceleration 
a' = a/(a^)^/^. The symbols are data from the LPT ex- 
periment at R\ = 460: x and + are the two measurements of 
the radial component and o are the axial component of accel- 
eration. The solid lines are the previous measurement of the 
radial component of acceleration using silicon-strip detectors 
in the same apparatus at R\ = 690 [2^. The inset shows 
a+*P(a+). 



for DLL{r), Eq. Q is solved for n(r) using Green's func- 
tions. The acceleration correlations i?LL(r) and RNNir) 
are then obtained from Eqs. pjl. 

We carried out three-dimensional Lagrangian Particle 
Tracking experiments in a von Karman water flow be- 
tween counter rotating disks. Our experimental tech- 
nique and particle-tracking algorithm have been de- 
scribed in detail previously [l2l 27 . [28j . Here, we re- 
port measurements from four experiments with Taylor 
microscale Reynolds numbers ranging from 200 to 690. 
The relevant parameters of the flow and the experiments 
are shown in Table [H All measurements were done in 
the same apparatus described in Ref. [5j, except for the 
R\ = 460 experiment, which was carried out in a new 
apparatus with a similar geometry but a different disk 
propeller. As a result, the integral length scale L of 
this experiment is markedly different from L in the other 
three. In the R\ — 460 experiment, we used the Phan- 
tom v7.3 cameras from Vision Research Inc., which are 
capable of recording at 37,000 frames per second at a res- 
olution of 256x256 pixels, nearly a 40% increase in frame 
rate compared to the v7.2 cameras used in the other ex- 
periments. Therefore, the R\ — 460 experiment has the 
highest temporal resolution among the four experiments 
reported here. 

We have shown before that the probability density 
functions (PDFs) of acceleration measured in the R\ = 
690 experiment agree well with previous measurements 
using silicon-strip detectors [13] ■ Figure [1] compares the 
acceleration PDF measured in the R\ = 460 experiment 
with the PDF measured in Ref. ^ a,t Rx = 690 using 
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43 
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690 


0.42 


1.2 


62 


30 


0.90 


24 
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TABLE I: Parameters of the experiments, u' is the root-mean-square velocity, e is the turbulent energy dissipation rate per 
unit mass. L = u'^/e is the integral length scale, rj and are the Kolmogorov length and time scales, respectively. Nf is the 
frame rate of the cameras, in frames per r^. The measurement volume is nearly a cube in the center of the tank, and its lateral 
size is given in the units of the Kolmogorov length scale rj. Ax is the spatial discretization of the recording system. The spatial 
uncertainty of the position measurements is roughly O.lAa;. Ns is the total number of acceleration measurements. 



silicon-strip detectors. (We note that there is possibly a 
weak dependence of acceleration PDF on Reynolds num- 
ber [5l.[3g|. Previous experiments, however, indicate that 
the dependence is so small that within experimental un- 
certainty, the measured acceleration flatness is nearly a 
constant over the range 285 < R\ < 970.) It can be 
seen that with a temporal resolution comparable to the 
silicon-strip detector measurement, where the sampling 
frequency is 65 frames per r,j, the PDFs measured with 
cameras are in remarkable agreement with the silicon- 
strip detector data. Even the fourth moment agrees 
well. The spatial resolution in the measurement with 
cameras, 40 /.tm/pixel, is significantly worse than that 
of the silicon-strip detector measurement (8 /im/pixcl). 
The agreement between the two measurements, however, 
suggests that using multiple cameras to determine the 
3D particle position results in better accuracy than the 
one-dimensional silicon-strip detector measurements. 

In Figure [21 we compare the longitudinal and trans- 
verse acceleration correlation coefficients measured at dif- 
ferent Reynolds numbers with the theoretical predictions 
by Obukhov & Yaglom 16], as obtained from Eqs. ([3]). 
It can be seen from Fig. 2(a) that the predictions agree 



well with the experimental data. As the Reynolds num- 
ber increases, the agreement between the predictions and 
the measurements increases. This is not surprising given 
that the viscous contribution to acceleration is neglected 
in the theoretical predictions. 

At large separations, the predicted acceleration corre- 
lations approach simple asymptotic scaling laws [l6j : 



RNN{r) 



^'-^2 ^3/2 ,-1/2 



-2/3 



(7) 

(8) 



Figure [2 (b) | compares Eqs. (O and ([8]) with experimental 
data. There are small but appreciable discrepancies be- 
tween the prediction and the measurements, which may 
be caused by the finite measurement volume and/or may 
refiect the need for still larger separations to see the 
asymptotic behavior. Another possible reason for the dis- 
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FIG. 2; (Color onhne) Comparison of the measured accelera- 
tion correlation functions with the Obukhov- Yaglom predic- 
tion [16] for (a) small separations and (b) large separations. 
The dashed and solid lines are the predicted longitudinal 
and transverse correlation functions, respectively. The sym- 
bols correspond to measurements from LPT data at different 
Reynolds numbers. Filled symbols are Rll and open symbols 
are Rnn. A - Rx = 200, V - Rx = 350, o - = 460, and 
□ - 7?A = 690. 
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FIG. 3: (Color online) Comparison of the measured pressure 
structure functions with the Obukhov-Yaglom prediction [l^ . 
The solid line is the prediction, and the dashed and dash- 
dotted lines indicate the r'^^^ and r'^^^ power-laws, respec- 
tively. The symbols are measurements from LPT data at 
different Reynolds numbers: A - Rx = 200, V - Rx = 350, o 
- i?A = 460, and U - Rx = 690. 



crepancy is that our flows are not isotropic. We observe 
anisotropy in acceleration even at the largest Reynolds 
number investigated, although the anisotropy decreases 
with increasing Reynolds number |5j| . 

Recently, Hill [3l[ proposed a refined theory for the 
acceleration correlations at small separations where the 
contribution from viscous forces is not neglected. Due to 
the limited spatial resolution in our experiment, however, 
our measurements cannot be used to test that theory. 

We obtain the pressure structure function by numer- 
ically integrating the equation for Rnn{t) in ® with 
experimentally measured transverse acceleration correla- 
tions. As already mentioned before, this equation holds 
in homogeneous, isotropic turbulence at high Reynolds 
numbers where the viscous contribution vanishes. 

In Fig. [3l we compare the prediction by Obukhov and 
Yaglom [Eqs. (j4]) to ([6])] with measurements obtained 
from the acceleration correlations. We plot the normal- 
ized pressure structure function Dp (r) = Ii{r) / p^ve. The 
measured data and the predictions are in good agree- 
ment over the range of separations accessible in the ex- 
periments. K41 inertial range scaling can be obtained 
from the limiting case of r 3> ?7 in the prediction of 11 (r), 
yielding Dp{r) ~ (r/r/)'*/'^, which is also plotted in Fig. [31 
This scaling law is close to the experimental data in the 
inertial range. The small deviations may be an indica- 
tion of intermittency. The finite Reynolds numbers and 
finite measurement volume effects could also contribute 
to the deviation in the scaling exponent. 

Recently, Bee et al. reported a r^/^ inertial range 
scaling for the pressure structure functions from DNS 
data up to Rx = 185. They also postulated that the r^/'^ 



scaling would persist at low Reynolds numbers and the 
K41 scaling might be observed only when R\ > 600. The 
r^/^ scaling law is shown in Fig. [3] for reference. As can 
be seen, all experimental data are much closer to the K41 
scaling rather than r^^^ scaling, and there is no appre- 
ciable change of slope over the range of 200 < Rx < 690. 
The discrepancies between experimental observation and 
numerical simulation remain to be investigated. 

Finally, let us note that the extent of the nominal iner- 
tial range measured from the spectrum or from structure 
functions can be very different, as has been shown for 
the case of Lagrangian velocity [12] . This subtlety could 
account for the difference between our experimental re- 
sults and previous investigations, in which the pressure 
spectra were studied. 

In summary, we simultaneously followed the trajecto- 
ries of multiple passive tracer particles in turbulent water 
flows with Reynolds number in the range 200 < Rx < 
690. The accuracy of the accelerations measured from 
the trajectories is comparable to previous single-particle 
measurements. We obtained spatial acceleration correla- 
tions from the multi-particle measurement and used the 
measured acceleration correlations to compute pressure 
structure functions from a relation that holds at high 
Reynolds numbers. We compared the measurements with 
theoretical predictions by Obukhov & Yaglom ^16^ and 
found that the predictions of both the acceleration corre- 
lations and the pressure structure functions arc in good 
agreement with the experimental data. We also observed 
K41 inertial range scaling in the measured pressure struc- 
ture functions over the range of Reynolds numbers inves- 
tigated. 
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